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We consider the equation of state of QCD at high density and zero temperature in perturbation 
theory to first order in the coupling constant a s . We compute the thermodynamic potential including 
the effect of a non-vanishing mass for the strange quark and show that corrections are sizable. 
Renormalization group running of the coupling and the strange quark mass plays a crucial role. 
The structure of quark stars is dramatically modified. 



The investigation of the equation of state for cold and 
dense strongly interacting matter and its consequences 
for the possible phases of quantum chromodynamics 
(QCD) is still in its infancy. From the observational 
point of view, compact stars provide the most promising 
"laboratory" , since their central densities might be high 
enough to allow for the presence of deconfined quark mat- 
ter. These dense objects are usually called quark stars 
or hybrid stars 0, and some observables, such as their 
mass-radius relation, may be used to constrain the equa- 
tion of state for strong interactions. 

In this Rapid Communication we compute the thermo- 
dynamic potential for cold quark matter with two light 
(massless) flavors, corresponding to the up and down 
quarks, and one massive flavor, corresponding to the 
strange quark, in perturbation theory to first order in 
a s in the Ms scheme. In this fashion, we can easily in- 
clude modern renormalization group running effects for 
both the coupling and the strange quark mass. We find 
that the corrections due to the running nonzero mass 
are sizable, and should not be neglected in the evalua- 
tion of thermodynamic quantities. Solving the Tolman- 
Oppenheimer-Volkov (TOV) equations using our equa- 
tion of state in the presence of electrons, we show that 
a running strange quark mass dramatically modifies the 
mass-radius diagram for quark stars, even at first order 
in a s . Although the numbers we present for the mass and 
radius of quark stars will be modified in a more realistic 
description, our results are consistent with constraints 
on the mass-to-radius ratio coming from measurements 
of gravitationally redshifted absorption lines in the X-ray 
burst spectra of the low-mass binary EXO0748 — 676 |2| • 

The thermodynamic potential for cold QCD in per- 
turbation theory to ~ a 2 s was first computed a long time 
ago N j 4 l |5 | . Nevertheless, the original approach to quark 
stars [7J, |S| made use of the bag model with corrections 
~ a s from perturbative QCD to compute the thermo- 
dynamic potential 22] . In the massless case, first-order 
corrections cancel out in the equation of state, so that 
one ends up with a free gas of quarks modified only by 
a bag constant. Finite quark mass effects were then es- 
timated to modify the equation of state by less than 5% 
and were essentially ignored. 



A few years ago, corrections ~ a 2 s with a modern def- 
inition of the running coupling constant were used to 
model the non-ideality in the equation of state for cold, 
dense QCD with massless quarks A perturbative 

calculation of the thermodynamic potential necessarily 
produces an unknown scale, A, associated with the sub- 
traction point for renormalization. In the case of zero 
temperature and mass and nonzero chemical potential, 
A is proportional to the quark chemical potential, \x. 
Nevertheless, on very general grounds, one can argue 
that reasonable values for A//i lie between 2 and 3, if 
one takes perturbative QCD as a model for the equation 
of state [lfj for cold strongly interacting matter. The 
former value corresponds to a strong first-order chiral 
phase transition, and the latter to a weak transition or 
a crossover. In the case of a strong first-order transition 
it was shown that there might be a new class of compact 
stars, besides the usual neutron stars, being smaller and 
much denser, with a large fraction of deconfined quark 
matter [1(1 lllj . This approach can be compared to treat- 
ments that resort to resummation methods and quasi- 
particle model descriptions 0, 0, Q, 0] . Remarkably, 
these different frameworks seem to agree reasonably well 
for ji >> 1 GeV, and point in the same direction even 
for [i ~ 1 GeV and smaller, where one is clearly pushing 
perturbative QCD far below its region of applicability. 

Even the most recent QCD approaches mentioned 
above generally neglected quark masses and the pres- 
ence of a color superconducting gap as compared to the 
typical scale for the chemical potential in the interior of 
compact stars, ~ 400 MeV and higher. However, it was 
recently argued that both effects should matter in the 
lower-density sector of the equation of state 0. In fact, 
although quarks are essentially massless in the core of 
quark stars, the mass of the strange quark runs up, and 
becomes comparable to the typical scale for the chemical 
potential, as one approaches the surface of the star. 

In what follows, we present an exploratory analysis of 
the effects of a finite mass for the strange quark on the 
equation of state for perturbative QCD at high density, 
leaving the inclusion of a color superconducting gap in 
this framework for future investigations. To illustrate 
the effects and study the modifications in the structure 



of quark stars, we focus on the simpler case of first-order 
corrections. Results including full corrections ~ a? s , as 
well as technical details of the calculation at each order 
and renormalization, will be presented in a longer publi- 
cation (lif . 

The leading-order piece of the thermodynamicpoten- 
tial of QCD for one massive flavor is given by 0, 0, 
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where N c is the number of colors and u = y /i 2 — m 2 . 

The first-order correction to Eq. Q, also known as 
the exchange energy 0, 0, , has the form 
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We work in the Feynman gauge and have adopted the 
abbreviation 
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and S(P — Q)= f35 Po , qo (2Tr) 3 5 3 (p — q), with temperature 
T = 1/(3. Here p = (2n + ljinT + [i is the frequency 
of the quark having four-momentum P and a chemical 
potential /z, and the sum is understood to be carried out 
over the integers n [6j. 

Using standard quantum field theoretical methods, one 
obtains the complete renormalized exchange energy for a 
massive quark in the Ms scheme (in the limit T — > 0): 
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The thermodynamic potential to order a s for one mas- 
sive flavor, given by the sum of Eqs. (@J and J[J, depends 
on the quark chemical potential fi and on the renormal- 
ization subtraction point A both explicitly and implicitly 
through the scale dependence of the strong coupling con- 
stant a s (A) and the mass m(A). The scale dependencies 
of both a s and m, which in the following we will take 
to be the mass of the strange quark, are known up to 
4-loop order in the Ms scheme |20j. Since we have only 
determined the free energy to first order in a s , we choose 
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FIG. 1: Q/Q,f at A = 2/i for fixed coupling (dash-dotted 
lines), including 2-loop running coupling (dashed lines) and 
including running of m 3 (full line). The light gray lines indi- 
cate m 3 = 0. 



where L = 21n(A/A M g), f3 Q = 11 - 2AT / /3, and (3 X = 
51 — IQNj/S and we take Nj = 3. The scale Aj-jg and 
the invariant mass rfi s are fixed by requiring \'2\\ a s ~ 0.3 
and m s ~ 100 MeV at A = 2 GeV; one obtains Aj-jg ~ 
380 MeV and m s ~ 262 MeV. With these conventions, 
the only freedom left is the choice of A. 

Note that Eq.QJ corresponds to an earlier result 01 
which was derived in a different renormalization scheme. 
Sticking consistently to the Ms scheme offers however 
some advantages, e.g. one choice of A allows one to use 
modern-day restrictions |2ll | on both the (/^-dependent) 
values of the strange quark mass and the strong coupling 
a s , whereas the authors of had to resort to (physi- 
cally motivated) choices for these. Also, our calculation 
is straightforward to extend to higher orders, since the 
coefficients for the renormalization group equations are 
known up to 4-loop order in the Ms scheme |2fJ. 

In Fig. ^ w e show the thermodynamic potential for 
one massive flavor (the strange quark), normalized by 
the free (massless) Fermi result, computed in the cases 
of fixed and running coupling and strange quark mass 
for A = 2/i. As one can see, renormalization group im- 
provements bring sizeable corrections in the region rele- 
vant for the structure of compact stars, \x < 1 GeV. At 
/i = 600 MeV, for which a s ~ 0.43, including a fixed 
mass for the strange quark decreases the thermodynamic 
potential by ~ 15%. The running of the mass (which 
insures renormalization scale invariance of Q up to order 
a s ) brings plus ~ 8% of decrease. 

To study the effect of the finite strange quark mass on 
the equation of state for electrically neutral quark mat- 
ter with 2 light (massless) flavors (up and down quarks) 
and one massive flavor (strange quark), we have to in- 
clude electrons, with chemical potential fi e . Imposing 



beta equilibrium, one has 
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Since neutrinos arc lost rather quickly, one may set their 
chemical potential to zero [lj, so that chemical equilib- 
rium yields 
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with (tx u , /id and Ms the up, down and strange quark chem- 
ical potentials, respectively. On the other hand, overall 
charge neutrality requires 
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where i%i is the number density of the particle species i. 
Together, the above equations insure that there is only 
one independent chemical potential, which we take to 
be m- Number densities are determined from the ther- 
modynamic potential by = — (dQ/dfi) and the to- 
tal energy density is given by e = fl + J2i^i n i^ where 
fl = YV + ^; 1 ' > ) arL d again i refers to the particle 

species. The pressure P is 
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where tlb — \ (n u + rid + n s ) is the baryon number 
density and the Gibbs potential per particle is given 
by = (/i«+Md + /-0- We restrict the freedom of 
choice for A(/x«, Md) Ms) by requiring that in case of van- 
ishing strange quark mass all quark chemical potentials 
and number densities become equal so that p( m ==°) = 
_Q(m a =0) anc ^ consequently, one has \x e — > JIs]- Fur- 
thermore, in order to compare our findings to exist- 
ing results in the literature 0, 0, we require that 
in the massless case A = 2m- Consequently, we choose 
A = | (m« + Mrf + Ms ) i but have tested that our results 
are not much affected by other choices obeying the above 
conditions 0. 

The effects of the finite strange quark mass on the total 
pressure and energy density for electrically neutral quark 
matter (plus electrons) are given in Fig. [21 There we 
show results for 3 light flavors and running coupling, cor- 
responding to the case considered in 0, and for 2 light 
flavors and one massive flavor, with both running cou- 
pling and strange quark mass (which reaches m s ~ 137 
MeV at m = 500 MeV). 

As can be seen from this Figure, there is a sizable dif- 
ference between zero and finite strange quark mass pres- 
sure and energy density for the values of the chemical 
potential in the region that is relevant for the physics 
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FIG. 2: Pressure and energy density scaled by Fermi values 
for A = | (/i u + fid + Ms)- We show results without renor- 
malization group improvements (dash-dotted lines), running 
coupling (dashed lines), both for m s = 0, and results with 
running mass and coupling (full lines). 



of compact stars, as was already expected from the re- 
sults presented in Fig. ^ As has been noticed by sev- 
eral authors 0,0,0], the resulting equation of state, 
e = e(-P), can be approximated by a non- ideal bag model 
form 



e = AB ef f + aP 



(10) 



Here a ~ 3 is a dimcnsionless coefficient while B e f f is the 
effective bag constant of the vacuum. Concentrating on 
the low-density part of the equation of state, one finds for 
massless strange quarks the parameters -B^j ~ 117 MeV 
and a ~ 2.81 while the inclusion of the running mass 
raises these values to -B^// — 137 MeV and a ~ 3.17 (all 
values having been obtained by including a running a s in 
the equations of state). Therefore, we expect important 
consequences in the mass-radius relation of quark stars 
due to the inclusion of a finite mass for the strange quark. 

The structure of a quark star is determined by the so- 
lution of the TOV equations Q] . One solves TOV by inte- 
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FIG. 3: Mass-radius relation for quark stars based on an equa- 
tion of state with running mass and coupling (full line) and 
for massless strange quarks (dashed line), for running cou- 
pling and with A = § (p u + + Mi- 



grating from some central starting pressure up to the sur- 
face of the star, where the pressure vanishes. The numer- 
ical solution of these equations in the case of our equa- 
tion of state for electrically neutral quark matter gives 
the mass-radius relation shown in Fig. [3J One can see 
from the Figure that corrections to the mass and radius 
of quark stars due to a running strange quark mass can 
be very large, ~ 25%. Also, while the most massive star 
for the m s = equation of state (with M/Mq ~ 3.2 and 
radius ~ 17 km) has a central density of tib — 0.5 fill - , 
this number increases to ub — 0.83 fm -3 (at /i = 470 
MeV) for the heaviest star (M/Mq ~ 2.16 at ~ 12 km) of 
the massive equation of state. One can put constraints on 
the mass-to-radius ratio by measuring the gravitational 
redshift at the surface of a neutron star. Our results 
for the mass-radius relation arc consistent with measure- 
ments reported by Cottam et al. 2], who obtain a red- 
shift z = 0.35 for the binary EXO0748 - 676, implying 
a ratio (M/Mq)/ R ~ 0.15/km, also shown in Fig. 
The inclusion of ~ a 2 , corrections to the pressure will in- 
crease its non-ideality and produ ce q uark stars which are 
smaller, denser and less massive . 

We have calculated the thermodynamic potential for 
cold quark matter with two light flavors and one massive 
flavor in perturbation theory to first order in a s in the 
Ms scheme. Results including full corrections ~ a 2 will 
appear soon 0] and will certainly modify the numbers 
we present here, since we know that corrections are large 
for this region of chemical potential [ToL Il3| . Moreover, 
for a realistic description of the structure of quark stars, 
one has to match the high-density equation of state onto 
a low-density hadronic one. This will also shift the values 
of masses and radii for the stars and, in some cases, allow 
for a new class of compact objects. Contributions due to 
color superconductivity |l6| as well as chiral condensa- 



tion 0| also affect this picture . Finally, if one aims 
for quantitatively tenable results, also the ~ a 3 contri- 
butions to the pressure should eventually be calculated. 
Therefore, our numbers for the maximum mass and ra- 
dius are meant only to be suggestive, and the point that 
we want to stress is the size of the corrections coming 
from the inclusion of a finite strange quark mass. 

In this Rapid Communication we present a still rather 
incomplete picture of cold and dense quark matter and 
the role it plays in the core of compact stars. Even so, we 
believe that our main conclusion is qualitatively correct: 
modifications in the pressure due to the inclusion of a 
nonzero strange quark mass are sizable and can dramat- 
ically modify the structure of quark stars. 
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